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(N Abstract 

We compute the lattice spacing corrections to the spectral density of the Hermitean Wilson 
Dirac operator using Wilson Chiral Perturbation Theory at NLO. We consider a regime where the 
quark mass m and the lattice spacing a obey the relative power counting m ~ aAq CD : in this 
situation discretisation effects can be treated as perturbation of the continuum behaviour. While 
this framework fails to describe lattice spectral density close to the threshold, it allows nevertheless 
to investigate important properties of the spectrum of the Wilson Dirac operator. We discuss the 
range of validity of our results and the possible implications in understanding the phase diagram 
of Wilson fermions. 
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I. INTRODUCTION 



Simulations with Wilson fermions at light quark masses are nowadays feasible thanks to 
considerable algorithmic improvements [IHS] developed in the last years. Since chiral sym- 
metry is explicitly broken, the spectrum of the Wilson Dirac operator is not protected from 
arbitrarily small eigenvalues, which might induce instabilities in numerical simulations [6]. It 
is therefore very important to have a theoretical understanding of the properties of the low- 
end spectrum of the Wilson Dirac operator. Moreover, spectral observables can be efficiently 
used to extract relevant quantities such as the quark condensate, like recently implemented 
in 0, yielding a further strong motivation to investigate the impact of lattice artefacts on 
the eigenvalues spectrum. 

In the continuum and in an infinite volume, the (renormalised) spectral density of the 
massive Dirac operator has a threshold given by the (renormalised) quark mass. Lattice 
artefacts are expected to change both the location of the threshold and the shape of the 
spectral density close to the threshold. Moreover, in a finite box of volume V one expects 
that when 'fHV ~ 1, where 7 is an eigenvalue of the massless Dirac operator and £ the chiral 
condensate, finite-size effects become important and can also induce relevant deformations 
of the spectral density with respect to the infinite-volume case. 

From the numerical point of view, some information about the low-end spectrum is pro- 
vided by lattice simulations with Wilson fermions carried out in the past years. For Nf = 2 
unimproved Wilson fermions, empirical observations |6] indicate that the median of the 
spectral gap distribution is linearly proportional to the quark mass m, while the width is 
basically independent on m and scales like ~ a/y/V, where a is the lattice spacing. On the 
other hand, for the 0(a)-improved theory, the situation is less clear [8] and those properties 
have not been confirmed. In [7] the mode number of the 0(a)-improved Wilson Dirac op- 
erator is computed, finding a nearly linear behavior up to ~ 100 MeV above the threshold. 
At the low end of the spectrum, a significant deviation from the continuum expectation is 
observed. 

From the theoretical side, Wilson chiral perturbation theory [§1 HO] (W^PT) is the tool 
which provides a systematic description of low-energy properties of lattice Wilson QCD 
including the leading discretisation effects. When approaching the chiral limit at finite lattice 
spacing, one enters in the regime where m ~ a 2 Aq CD , which is where discretisation effects 
compete with the quark mass to the explicit breaking of chiral symmetry. Lattice artefacts 
induce a non-trivial phase diagram, and two different scenario have been foreseen: in the 
so-called Aoki scenario there is a range of quark masses (in the Aoki phase [TT] ) where there 
are two massless pions. On the other hand, in the so-called Sharpe-Singleton scenario [9] 
there is a first-order phase transition, and the three pions remain massive in the chiral limit. 
From the point of view of W^PT, in this regime the discretisation effects appear already at 
leading order in the chiral expansion. Lattice artefacts in the infinite volume spectrum of 
the Hermitean Wilson Dirac operator have been computed in this regime in [12J, although a 
working framework has been only found by assuming additional conditions on the couplings 
of W^PT associated to discretization effects. 
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In a recent study [131 E] the same power counting for a has been adopted, but in a 
finite-box in the e-regime 1 . In this case the LO predictions of W%PT can be obtained also 
by means of a "modified" Random Matrix Theory which includes 0(a 2 ) effects. The main 
difficulty of this computation is that it involves exact integrals over the zero modes which 
must be defined at fixed topological charge introduced via the number of real modes of the 
Dirac operator, which nevertheless has an intrinsic ambiguity for Wilson fermions at finite 
lattice spacing. 

In this work we will study the discretisation effects in the spectrum of the Dirac operator in 
a different regime. We consider W%PT with the power counting m ~ aAq CD (GSM regime), 
where lattice artefacts appear only at next-to-leading order in the chiral effective theory, 
and can hence be treated as a perturbation. This simplifies considerably the computation; 
it allows nevertheless to extract important information about the impact of lattice artefacts 
in the spectral density. 

This paper is structured as follows: in sec. [TT] we recall definitions and properties of the 
spectral density of the Wilson Dirac operator; in sec. Ill we present the setup for W%PT 
needed for the computation; in sec. IV we give the details about the calculation of the 
spectral density, including the finite- volume corrections (in the p-regime); in sec. [V] we 
discuss our results and show a comparison with numerical data. All technical details about 
the computation are deferred in the Appendices. 



II. SPECTRAL DENSITY FOR THE WILSON DIRAC OPERATOR 

We start considering QCD for Nf = 2 degenerate quarks in infinite volume with current 
quark mass m. The average spectral density of the massless Hermitian Dirac operator — iD 
can be defined as 

p D (7, m) = Jim — ^(5 (7 - Jk)}, (1) 

k 

where 7^ are the eigenvalues of the massless Hermitian Dirac operator and (...) indicates 
the usual path-integral average. The Banks-Casher relation [18] tells us that the spectral 
density is related to the chiral condensate S in the following way 

PD (7, m) = - [1 + 0({| 7 |, m}/A Q c D )} ■ (2) 

7T 

One useful way to determine this relation in continuum QCD [191 120] consists in adding a 
valence quark ip v of mass m v to the theory. Using the spectral decomposition for the valence 
chiral condensate one obtains 

®M = -hi p -^. (3) 



1 Standard W%PT has been extended to the e-regime in [T5lU7| . 
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This relation can be inverted because the spectral density is independent on the valence 
quark mass 

Disc [{ip v ip v }] \m v =- iy = -2^(7, rn), (4) 

where Disc indicates the discontinuity across the imaginary valence quark mass. Strictly 
speaking eq. ^ is ultraviolet divergent also after renormalisation of the quark masses and 
the gauge coupling constant and one must introduce a cutoff in the integration range. The 
ultraviolet divergences turn out to cancel when computing the discontinuity of the valence 
scalar condensate in eq. Q. Thus it is enough to compute the valence condensate for real 
masses and analytically continuing the resulting expression for complex masses. Eq. Q 
can be used naturally in chiral perturbation theory and one obtains [20] the Banks-Casher 
relation ^ and the NLO corrections [71 12~T]. 

We now discretise our Nf = 2 continuum QCD action on a lattice of spacing a and we 
consider Wilson fermions. All our considerations can be generalised to the case of Wilson 
twisted mass. To study the spectral density of the Wilson operator one has to take into 
account that the Wilson operator D\y is neither Hermitian nor anti-Hermitian and it has 
complex eigenvalues. It is thus advantageous to define the Hermitian Wilson-Dirac operator 

m 

Qm = 75 {D w + m) = Ql. (5) 

We indicate now the spectral density of Q m as Pq(X, to,) and the one of Q 2 m as p(a, to), where 
A and a label the eigenvalues of the corresponding operators. 

If the continuum theory has been regulated on a lattice the ultraviolet divergences of the 
chiral condensate appear as power law and logarithmic divergences in the lattice spacing a. 
Moreover, when using Wilson fermions the lack of chiral symmetry implies that the scalar 
condensate is not directly related to the spectral density as in the continuum. It is thus not 
obvious how to work out the renormalisation and 0(a) improvement of the spectral density 
starting from eq. Q or from the corresponding version for Wilson lattice QCD. To study 
the renormalisation and the 0(a) improvement of the spectral density one needs to relate 
it to correlation functions of local operators where standard arguments on renormalisability 
and 0(a) improvement can be applied. 

It has been shown in ref. [7J, using chain correlators of scalar and pseudoscalar densities 
how the spectral density of Q 2 m renormalises. In the following we will be mostly interested 
in the spectral density of Q m that renormalises as follows 

[Pq] r (X, m R ) = Z P p Q (ZpX, to), (6) 

where rriR is the renormalised quark mass. Moreover, additionally to the standard improve- 
ment of the action and local operators there are 0(am) cutoff effects which need to be 
removed to fully improve the spectral density [7J . This analysis guarantees that with Wilson 
fermions the spectral density is a well defined renormalisable quantity and with a well de- 
fined Symanzik expansion. We can thus compute the spectral density in the chiral effective 
theory using a generalisation of eq. Q for Wilson chiral perturbation theory. 
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A further spectral observable which can be defined and measured in a lattice simulation is 
the integrated spectral density 

N(A 1 ,A 2 ,m) = dX p Q (X,m), A 2 > Ai > m, (7) 

JAi 

which represents the density of modes in the interval between Ai and A 2 ; it contains the same 
physical information as the spectral density and satisfies N R (A 1R , A 2) r, m R ) = N(Ai, A 2 , m), 
i.e. it is a renormalisation-group invariant [TJ. 

To relate a partially quenched condensate with the spectral density of Q m it is convenient 
to introduce a doublet of degenerate twisted mass fermions Xv with a mass term ip v j 5 T 3 , 
where r 3 is the third Pauli matrix. The relation between the condensate and the spectral 
density is now different but can be worked out and it reads |12j 

&7srV> = f JA W(*."> ( 8 ) 



oo 



The need to symmetrise in A the spectral density in the numerator of (|8j) can be understood 
from the lack of chiral symmetry of the Wilson operator which renders the spectral density 
not symmetric when A — > — X. Inverting eq. (Isb one obtains 



Disc [{Xvls^Xv}] L=a = 2iir [p Q (X, m) + p Q (-X, m)] (9) 

where the untwisted mass in the valence sector coincides with the sea quark mass m. Eq. Q 
is our starting point for the computation of the spectral density in Wilson chiral perturbation 
theory. 

We stress that in chiral perturbation theory the problem of power law divergences is 
absent because those are mapped to the presence of the so called high-energy constants 



usually indicated with H. As we will see in sect. [TV] this is indeed the case for the expectation 
value of the pseudoscalar condensate in the chiral effective theory. This implies that in 
the effective theory the computation of the discontinuity is a well defined procedure. It 
remains to make sure that the spectral density computed in the effective theory through the 
discontinuity of the partial quenched chiral condensate does not overlook additional cutoff 
effects which go beyond the standard 0(a) corrections stemming from the action or the 
local operators. Our previous discussion and the results of ref. [7] guarantees that the only 



additional cutoff effects are of Of am) and as we will see in sect. Ill these cutoff effects are 



of subleading order in our power counting of Wilson chiral perturbation theory. 



III. SPECTRAL DENSITY IN WILSON CHIRAL PERTURBATION THEORY 

In order to investigate the discretisation effects in the spectral density of the Hermitian 
Wilson Dirac operator we compute the valence pseudoscalar density in the framework of 
Partially Quenched Wilson chiral perturbation theory (PQW^PT). According to eq. (J8|, we 
consider Nj = 2 sea quarks with bare mass m and a doublet of valence twisted mass quarks 
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with (m + ifi v r 3 ). We formulate the effective theory on a finite volume V = L 3 T and we 
keep the quark masses in the p-regime, corresponding to the power counting 

m,/j v ~0(p 2 ), 1/L,l/T~0(p) (fO) 

in terms of the momenta p. The extension of the effective theory to the case of non-zero lat- 
tice spacing is done in two steps: after matching the lattice QCD action with the appropriate 
Symanzik continuum effective action [221 [23] . one writes down a chiral Lagrangian which 
contains the standard continuum terms plus additional operators that transform under chi- 
ral symmetry as the operators of the Symanzik effective theory. For the Wilson action, this 
has been studied in [9], [10] . When introducing lattice artefacts in the chiral effective theory, 
we have to define as usual the relative power counting between the quark mass and the 
lattice spacing a. In this work we adopt a counting corresponding to the so-called GSM 
(Generically Small quark Mass) regime [24"l [25] . where m, \l v ~ aAq CD : in this case the 
explicit breaking of chiral symmetry is dominated by the quark mass, and lattice artefacts 
can be treated as perturbations. 

Partially quenching can be implemented in the Chiral Effective theory by means of two 
techniques, namely the graded-symmetry method [26,, 27] and the replica method [28] . In the 
first one, one introduces "ghost" quarks, whose determinant cancels the one of the valence 
quarks. In the second one, one enlarges the valence sector to N r flavors and one eliminates 
the corresponding determinant by sending N r — > 0. The equivalence of the two methods has 
been shown at the perturbative level in [28] ■ In this section we will focus on the graded- 
symmetry method; we checked nevertheless that the same results can be obtained with the 
replica method, which we summarise in App. [B} For our specific case we have to consider a 
chiral Effective Theory with a graded symmetry group SU(4|2) L xSU(4|2)# spontaneously 
broken to SU(4|2)ij + £. The formulation we adopt for our computation using the SU(4|2) 
effective theory is based on an extension of the framework proposed in |29j for the SU(3|1) 
case. We consider a mass matrix of the form 

Ai = diag( ^jn^ ,m + i/j v r 3 ,m + ifi v r 3 ). (11) 

2x2 sea 2x2 val 2x2 ghost 

The quark mass in the ghost sector has been already set equal to the mass in the valence 
sector. Moreover, the untwisted part of the valence quark mass has to be equal to the sea 
quark mass. We parametrise the pseudo Nambu-Goldstone bosons by the field U G SU(4|2) 

U(x) = u v e 2i ^ F u v , £ = ( 12 ) 

a 

where F is as usual the pseudoscalar decay constant, and T a , a = 1, . . . , 35 represent the 
generators of the corresponding Lie algebra, which satisfy 

ab 

Str(T a T fe ) = y —. (13) 



We refer to App. [A]for the explicit form of g ab (eq. ( A15 )), and for a summary of conventions 



and properties of the SU(m|n) group. The constant field uy represents the ground state of 



6 



the theory, which can be obtained by minimising the potential in the LO Chiral Lagrangian. 
The solution of the potential minimisation yields 2 

u v = diag( ^ , ^WJ, ^) , (14) 

2x2 sea 2x2 val 2x2 ghost 

where 

sina; = coswo = , m P = \ffn 2 + u?.. (15) 

m P m P 

In the language of twisted mass Wilson theory, mp is the so-called polar mass. 

Taking into account our power counting, the full NLO Chiral Lagrangian in PQW^PT 
can be written as 

C = C 2 + £ 4 + C a . (16) 

In the GSM regime, the LO Lagrangian can be written in the same form as the continuum 
one 

F 2 

C 2 = —{8tr(d fl Ud ll rt)-2BStT{Mrt + UMi)}, (17) 
provided we substitute the quark mass with the so-called shifted mass [9], which incorporates 



the leading 0(a) corrections. In the following we assume that m in eq. (11 ) includes already 
this shift. The coupling B is related to the chiral condensate E, B = S/F 2 . 
From £ 2 we can extract the propagator 

(C(x)e(y)) = 9 ah G\(x-y- M 2 a )+2B(m P -m)h ab G 2 v (x-y, M 2 ), (a, b = 1, ... , 35), (18) 
with 

GUx,M 2 ) = -y—^——, (r>l), p = 27r(^,^A^y n 1234 eZ. (19) 



We recall that [30] 

G r v (0,M 2 ) = G r (0,M 2 )+g r (M 2 ), (20) 

where 

^•''^w/^PTn' (21) 

is the infinite- volume contribution, while g r (M 2 ) represents the finite- volume correction and 
is UV-finite. A specific representation of g r (M 2 ) will be considered in App. [c] The matrix 



h ab is defined in App. |Aj eq. (A17). The mass term M 2 appearing in the propagator is given 
by 

r M 2 S = 2mB, a = 1,...,3 

M 2 = i M 2 V = (m + m P )B, a = 4, . . . , 11, 16, . . . , 23 (22) 
[ M 2 V = 2m P B, a = 12, . . . , 15, 24, ... , 35. 

2 In the graded symmetry formulation of PQ^PT there are subtleties concerning the minimisation of the 



potential. We will discuss them in more detail in sec. IV 
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The continuum NLO PQ Lagrangian reads [31] 

£ 4 = 2BL 4 Str(d fl U j d tl U)Sti(MU j + M j U) + 2BL 5 Sti(d fl U j d II U(MU j + M j U)) 

-4B 2 L e Str(U j M + M j U)Str(U j M + M ] U) (23) 
-4B 2 L 7 Sti{M*U - MU j )Str(M j U - MU ] ) 
-4B 2 L 8 Str{MU j MU ] + M ] UM ] U) - 4B 2 H 2 Str(M j M), 

where we have written down only terms which can enter in our specific computation. The PQ 
Chiral Lagrangian encoding discretisation effects up to 0(a 2 ) can be written as [TOl |32| [33] 

C a = -2BaW 6 Str{MU j + M j U)Sti{U + U j ) - 2BaW 7 Str{M j U - MU j )Str{U - U r ) 
-2BaW 8 Str{MU j2 + M ] U 2 ) - a 2 W^Sti(U + U ] )f - a 2 W 7 {Str{U - U ] )f 
-a 2 W£ti{U 2 + U ]2 ), (24) 

where a = 2W a has dimension [energy] 2 and W is the LO low-energy coupling absorbed 
in the shifted mass m. Also in this case we have disregarded terms which are not relevant 
for our computation. 

IV. CALCULATION OF THE SPECTRAL DENSITY WITH SU(4|2) GRADED 
GROUP METHOD 

In this section we compute the spectral density at NLO in W^PT, making use of the 
definitions and assumptions introduced in the previous section. The observable from which 
we start is the partially quenched twisted pseudoscalar condensate made of valence quarks 
defined in eq. (j8J). In order to extract it we have to introduce source terms in the Chiral 
Lagrangian by means of 

M M p , = M+ P3 f 3 , M*^M\ )3 =M* -p 3 T 3 , (25) 

where M. is the mass matrix introduced in eq. (11), and f 3 has non-zero components only 
in the valence sector 

f 3 = diag( , r 3 ^ , ). (26) 

2x2 sea 2x2 val 2x2 ghost 

The pseudoscalar density is obtained by deriving the action associated to the Chiral La- 



grangian in eq. (16) with respect to p% 



v3(x) = ^y 5| -»- (27) 

At LO (0(p°)) we obtain the (continuum) expectation value 

(V 3 ) LO = 2iSsinw - (28) 
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The NLO result is 



(V 3 ) NLO 



2iE sin uq < 1 + 5 cot coq + 



F 2 



-G l v (0, M 2 J - 2(4(0, M 2 SV ) + 32L 6 Ml (29) 



- -(M£ v - M 2 s )G v (0, M,fJ + 4M^(tf 2 + 2L 8 ) + 8a (2W 6 + W 8 cosu ) 



where 5 is an 0(p 2 ) quantity that represents the correction to the ground state angle uq due 
to lattice artefacts. In the following subsection we discuss the computation of this effect. 



A. NLO correction to the vacuum 

The 0(p 2 ) terms in the chiral Lagrangian give rise to a shift in the vacuum angle uq [25] . 
which must be computed by minimising the NLO potential. An important fact is that at this 
order the continuum Lj do not contribute to this realignment, and only the discretisation 
effects must be taken into account. Since this is an important point of the calculation, we 



start by recalling how the continuum ground state given in eq. ( 14 ) is obtained. Using the 
equations of motions one can show that uy has to commute with the matrix Ai^Ai, which 
implies that uy is diagonal. Therefore we can parametrise it as 

uy = diag ( e 2 , e 2 , e 2 , e 2 , e 2 , e 2 j ; (30) 



where 7r < 1; . . . , 4 < tc are standard phase factors. The condition Sdet(wy) = 1 requires 

01 + 02 + 03 + 04 - 05 - 06 = 0. (31) 

As suggested in [34] . we perform an analytic continuation for the ghost components 

05 = «05, 06 = ^06, (32) 

where now 5 6 are real variables taking values along the real axis. The LO potential can be 
separated in quark and ghost components 

V LO = V q , LO + V g , LO , (33) 

with 

V q ,LO = 2[mp(cos(w - 03) + cos(u;o + 04)) + m(cos0i + cos0 2 )], (34) 

V g ,LO = -2rnp[cos(u) — ife) + cos(u) + i<fi 6 )]. (35) 

The potential is complex, and a minimisation procedure in the usual sense is not applicable. 
The minimisation of the potential in Euclidean field theory is equivalent to perform a saddle- 
point expansion of the functional integral, and the saddle-point expansion can be performed 
even if the potential is complex. Our prescription, following [35], is to find the saddles for 
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complex 5 6 and then deform the contour of integration for each field variable in order to 
pass through the saddle points. The saddle points are chosen following two criteria: a) in 
order to maximise the value of the real part of the potential at the saddle; b) in order to 
have the direction where the real part of the potential rises more steeply consistent with the 
possibility of deforming the contour integration without encountering the other saddle points. 
This procedure has been discussed in detail in [35] in the context of the phase diagram of 
quenched Wilson ChPT. We find that both criteria for the choice of the appropriate saddle 
points are fulfilled by 05 = —06 = uq. In this way one obtains the expected result given in 
eq. @4}. 

We now repeat this procedure including the NLO lattice artefacts. We consider a ground 
state uy, nlo of the form given in eq. (30) and we perturb the LO solution by choosing 

0i = 02 = 5 S , (36) 
03 = + 5 V , 04 = -(wo + 5 v ), (37) 
05 = u + 5 g , 6 = -(u + 5 g ), (38) 

where S S)V>g are corrections of 0(p 2 ). We then perform the same analytic continuation on 
the ghost components, 

05 = ^05 = ^(^0 + 4), 06 = i(f>6 = -i{uo + S g ). (39) 

By expanding the NLO potential at 0(p 4 ) we obtain 

Vnlo = «o + ai$s + a 2 8s + a 3 5 v + a 4 5 2 v + a 5 8 g + a 6 5 2 g + . . . = V q%NLO + V g>NLO , (40) 

with 

a = -2 (2a 2 (AW^ + W' & ) + m{AaB{AW & + W 8 ) + £)) , (41) 
ai = 
a 2 = mS, 

a 3 = ia 5 = 16a sin w (aW$ cosw + 2aW$ + BmpW 8 + 2BmW e ), 
cla = ci6 = mpTi. 

Also in this case the potential can be split in a real (sea + valence) part V q ^LO and a 
complex (ghost) part V 9 ^lo- The saddle point expansion yields 

o s = 0, 6 v = -p-, 5 g = -^ = -i5 v . (42) 
2a 4 a 2a 6 

The ground state of the theory at NLO is then given by 

U V , NLO = ( ^ , £^fV2, ^_W2) , (43) 

2x2 sea 2x2 val 2x2 ghost 

with 

. . . 16asinw f /„. W 8 \ 2a (. W' % \\ 

5 = 5 v = 5 g = — 2 | (W 6 cos w + ^)+Jp[ W 6 + ^ L cosw ») ) • ( 44 ) 

Since the discretisation effects represent a perturbation in our power counting, this solution 
does not impose constrictions on the absolute value or the sign of the LECs W$, W 8 , W$, W$. 
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B. Spectral density in infinite volume 



The spectral density can be computed from (V 3 ) by taking the discontinuity along imag- 
inary p v , as given in eq. (J9]). At LO we obtain the known result 



[p Q (X,m)+ p Q (-X.m)] 



2£A 



LO 



71 



(45) 



which does not depend on the volume and on the lattice spacing in our specific regime. The 
NLO result in infinite volume is given by 



[p Q (X, m) + p Q {-\, m)] NLO 



2£A 



m 2 A 



7T 



+ 



A 2 -m 2 (4vr) 2 F 4 



7T 



v 7 ^ 



+ m(3LQ — 1) + 2V 7 A 2 — m 2 arctan 



ni 



(46) 



where A arises from the correction to the ground state and is given by 



~ 16a / XT 2aW' 

A = We, H 

F 2 V Ml 



If we define a shifted sea quark mass as follows 



(47) 



(48) 



under the condition that A ^> rh, the correction proportional to A can be resummed. In this 



rh = m ( 1 + A ) , 



case one could then rewrite eq. (46) by omitting the term proportional to A and substituting 



m — ?• rh everywhere. This resummed formula would be equivalent to the one in eq. (46) up 
to higher order corrections. 



Notice that the term ajM 2 ss in eq. (47) is not singular in the GSM regime. The UV diver- 
gences arising from G 1 (0, M 2 ) and G 2 (0,M 2 ) have been cancelled by defining renormalised 
couplings [36]. Using dimensional regularisation in 4 — 2e dimensions and adopting the 
convention of [7] we define 



3/i 



-2c 



64(4vr) 



L 6 h 7 — In An — 1 

e 



(49) 



where fi is a renormalisation scale. 

The integrated spectral density A^(Ai,A 2 ,m) defined in eq. ^ can then be computed 
straightforwardly. For better readability we report it in App. [Cj eq.( |Cl[ ). We stress that in 
the continuum these expressions coincides with the ones computed in [7]. 
It is useful to express our results for the spectral density in terms of the PCAC quark mass. 
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The NLO lattice corrections to the PCAC mass can be computed in the conventional W%PT 
for Nf = 2. The result is [25] 

16a W 8 W w 2a ( W£ 



m PCAC = m{l + — + ^ + ^ + W + ) h (50) 



where W\q is an extra LEC. Eq. (46) can be then rewritten by substituting everywhere 
m mpcAc and A — y A, with 

A + ^ + m 



F 2 \ 2 4 M 2 S/ 
As before, the correction given by A can be resummed; the shift in the quark mass given in 



eq. (48) can be rewritten as a shift in the PCAC mass: 

(l + A) = mp CAC (1 + A) . (52) 



m = m 



Notice that at NLO we can consistently substitute m — » m PC Ac m M% s 
In fig. [Owe plot [pq(X, m) + Pq{— A, m)]jvz,o in the continuum theory and in the discretised 
case with We = Ws = W\o = 0, corresponding to a non-perturbatively 0(a)-improved 
theory. We consider in particular a quark mass m PC Ac = 26.5 MeV; the values of other 
parameters used in this plot are specified in the caption and justified by our numerical 
analysis discussed in sect. |Vj 

The central black curve corresponds to the continuum case; W$ > (A < 0) gives rise to 
the red curve, while W$ < (A > 0) corresponds to the blue curve. The dashed lines for 
A < 40 MeV are a reminder that our result is not expected to be valid close to the threshold. 
This issue will be discussed in more details in sec. [V] 



C. Finite volume corrections 

Finite-volume effects arise at NLO and can be computed by taking into account the 



corrections to the propagators given in eq. (20). We can write down the spectral density as 



sum of the infinite volume and finite-volume corrections: 

Pq(X, m) + pq(-\, m) = p Q (\, m) + p Q (-X, m) + Ap^(A, m) + Apg(-A, m). (53) 
The explicit expression for the finite-volume correction Apg(A,m) + Ap^(— X,m) is given 



in App. |Cj eq. (C4). By integrating over A we obtain the corresponding finite- volume 
corrections for the integrated density 

N V (A U A 2 , m) = JV(Ai, A 2 , m) + AA^ y (A 1; A 2 , m). (54) 

In eq. (C5) we give the formula for AiV y (Ai, A 2 , m) for the particular case Ai = m. 



We stress that in our power counting finite-volume effects are insensitive to lattice arte- 
facts at NLO, and these results are the one of the continuum theory already obtained in [7]. 3 

3 The formulae for finite- volume effects are not given explicitly in [7| ; we checked nevertheless that our results 
are in agreement with them |29j and we decided to report them in the present work for completeness. 
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FIG. 1: The spectral density [pQ(X,m) + pq(— A, m)]NLO in the infinite volume. We used the 
parameters S = (275 MeV) 3 , m PC AC = 26.5 MeV,F = 90 MeV, L 6 = 5, p = 139.6 MeV. The solid 
black line corresponds to the continuum xPT prediction, while the red (blue) lines correspond to the 
lattice W^PT prediction (O(a)-improved) on eq. (46) with a 2 Wg = ±5 • 10 6 MeV 4 , corresponding 
to A = TO. 53. 



In fig. Swe show the relative correction AN v (m, A, m)NLo/N(m, A, m) nlo as a function 
of A, for two different volumes with T = 3.84 fm and L = 1.92 (red dashed curve) and 
L = 2.56 fm (black solid curve). 

Notice that the finite- volume correction N v (Ai, A2, to) diverges at the threshold (Ai = A2 = 
to): this could be interpreted as a signal that, even if the quark mass is in the p-regime, 
another scale comes into play in this problem, namely v A 2 — to 2 , and it needs to be treated 
with a different power counting when approaching the threshold at finite volume. Since in 
the e-expansion of the chiral theory divergences naturally appear when taking the chiral 
limit at finite volume, this might be the correct power counting to adopt for this additional 
scale when y/X 2 — m 2 YM ~ 1. Having said that, for L > 2.0 fm and sufficiently far from the 
threshold, finite-volume corrections amount to few percents. 



V. DISCUSSION OF THE RESULTS AND CONCLUSIONS 



The important results of this paper are the expression of the spectral density eq. (46) 
and its integrated form (CI). They are both obtained in the p-regime of W%PT at NLO. 
Eq. (46) shows that the spectral density, when computed on a lattice with Wilson fermions, 
is modified in two aspects with respect to the NLO continuum formula (see for example 
refs. [37] or simply set all the H^s in eq. (46) to zero). The lattice artefacts modify the 
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FIG. 2: The finite-size effects for the (continuum) integrated spectral density, 
AN V (m, A,m)NLo/N(m, A,m)NLo as a function of A. We used the parameters £ = (275 MeV) 3 , 
m = ttt-pcac = 26.5 MeV, F = 90 MeV, Lg = 5, n = 139.6 MeV. The black solid line corresponds 
to L = 2.56 fm, T = 3.84 fm, while the red dashed line corresponds to L = 1.92 fm, T = 3.84 fm. 



behaviour of the spectral density near the threshold, i.e. the term proportional to A in (46), 
and and its absolute normalisation, i.e. the term proportional to W% in (46). The absolute 



normalisation correction proportional to Wq vanishes if the theory is non-perturbatively 
0(a)-imp roved. If the theory is not improved, beside the term proportional to Wq, the 
correction to the normalisation contains in principle an additional O(a) term stemming 



from the renormalisation constant Zp. We prefer to omit in our final formula (46) this term 



because it strictly depends on the way Zp is computed and only in very few cases we can 
have a suitable representation in W;\PT [3B] • We stress nevertheless that when analysing the 



spectral density computed in the unimproved theory using eq. (46) the O(a) cutoff effects 
of Zp have to be taken into account in some way. 

More interestingly the cutoff effects modify, the A and m dependence of the continuum 



formula. These corrections, expressed in terms of the PCAC mass (see eqs. (50 52)), depends 



on the following LECs: W%, W\q and W%. The modifications induced by cutoff effects to 
the continuum formula are plotted in fig. [T] To simplify the discussion we have chosen a 
lattice setup where the theory has been non-perturbatively improved, thus we set all the 
O(a) LECs Wq = Wg = Ww = 0. With the black curve we plot the A dependence of the 
spectral density in the continuum as given by NLO ;\TT [3 EZ] at given values for m, S 
and F (see the caption of the plot). In the same plot we also show the modifications to 
the continuum formula induced by W^PT at NLO. The first observation we make is that 
the way the spectral density depends on A, close to the threshold, is rather different from 
the continuum formula. The second observation we make is that the W^PT formula has 
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a rather different behaviour depending on the sign of Wg, especially close to the threshold. 
If A > when A — > mpcAc the W%PT formula has a non-integrable singularity, while the 
continuum formula has an integrable singularity. If A < the spectral density goes to zero 
for a value of A > mpcAc- These behaviours are shown in fig. [T] This A dependence of 
the spectral density near the threshold is totally dominated by the lattice NLO corrections 
induced by the vacuum realignment, indicating that our expansion cannot be trusted for too 
small A. Our power counting breaks down near the threshold: more specifically, it is valid 
only when \/A 2 — m? ~ ciAq CD , i.e. when the scale a/A 2 — m 2 obeys a GSM power counting. 
By moving towards the threshold at finite lattice spacing one enters in the region where 
a/A 2 — m 2 ~ a 2 AQ CD (Aoki regime), and lattice spacing corrections can not be treated in a 
perturbative way as it is done in this work. In fig. [T] we draw dashed line where we assume 
our formula for the spectral density to break down. The minimal A ~ 40 MeV is just an 
estimate obtained fitting the available numerical data, as we will discuss futher on. It seems 
a reasonable estimate given the fact that for A ~ 40 MeV the relative lattice spacing NLO 
corrections are still sufficiently small. Of course this is just an estimate based on the value 
of A we input, i.e. from the size of the lattice artefacts. Only a thorough analysis of the 
numerical data can determine the value of A, thus the range of validity of our results. 



Thus we can try to use our formula (CI) to fit the available numerical data and see if 
we obtain reasonable estimates of £ and A. To compare our formula with the available 
numerical data from ref. [7j we decided to integrate the spectral density from A = Ai and 
A = A 2 , taking Ai = 40 MeV, in order to avoid the small A region where the lattice 
corrections to the xPT formula dominate. In order to compare with the results of [7J we 
use the mode number defined as z/(Ai,A2,m) = VN(Ai, A2, m) and we subtract the value 
v(0, Ai,m) from the numerical data. In fig. [3] we show a typical fit of the numerical data 
of ref. [7J using the infinite volume formulae of continuum xPT and W%PT in the range 
A2 G [40 : 60] MeV. The numerical data are obtained at a value of the renormalised PCAC 
mass mp, — 26 MeV. 

We have performed fits using the both the NLO formula? of %PT and W^PT. 4 With the 
continuum xPT formula we have fit £ and fixed the following parameters: fi = 139.6 MeV, 
F = 80, 90 MeV, I 6 = 3, 4, . . . , 7 and m PCA c = 26.5 MeV. With the W^PT formula we 
have set all the O(a) LECs W$ = W% = W^io — because the numerical data are non- 
perturbatively 0(a)-improved. Additionally we have fixed the value of X 1 / 3 within the 
range 250 — 300 MeV and we have fit A. We have performed fits with several forms for 
the W%PT formula, all differing by NNLO terms. We have observed that the best fits are 



obtained with eq. (CI). In order to decide on the 'best fits' we have selected the fit results 
that gives small values of the squared differences between the input data points and the 
theoretical formula evaluated at the same points. 

A summary of the fit results, satisfying our criterion is given in fig. |4| The plot shows 
the parameter space of the fit results in the A, E 1//3 plane. Each coloured band represents 
a fixed value of L 6 and the lower and upper borders of each band represent respectively 



A LO fit with %PT does not perform better in any range of A where the numerical data are available. 
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FIG. 3: Comparison of the numerical data (red curve) for the mode number z^(0, A 2 ,rn.) and the 
NLO fit results using continuum xPT (black curve) and W%PT (blue curve). The fits are obtained 
using the following parameters: Ai = 40 MeV fi = 139.6 MeV, F = 90 MeV, Lq = 5 and 
mpcAC = 26.5 MeV. The fit range is A 2 G [40 : 60] MeV. 



F = 80 MeV and F = 90 MeV. The two most important observations are that reasonable 
fits are obtained only for 4 < L 6 < 7 and for almost the whole parameter space we obtain 
a negative value of A. Our best fit result shown in fig. [3] gives A = —0.55 with a value of 
j]i/3 _ 273 MeV, which is in the right ballpark. We will come back to the consequences of 
this result at the end of this section. 

We remark that while the fit results obtained using the continuum NLO formula are not 
disastrous always a better agreement is found using our W%PT result. 
The W%PT formula is able to describe well the numerical data also in a larger A 2 range, i.e. 
[40 : 100] MeV, with reasonable values of the chiral condensate. 

Even if our analysis is rather qualitative and more numerical data are needed, we can try to 
draw some conclusions. W%PT describes the numerical data better than continuum xPT 
in the range of applicability of our results, i.e. sufficiently away from the threshold, but 
we cannot exclude that adding NNLO terms in the continuum formula would improve the 
quality of the fits. 

Our fit results seem to prefer a negative value for A implying a positive value for W%. In 
the literature, the LECs W$ and W$ are often found combined as [39] 



a 2 



^ 6 + T" J=-^6-' (55) 



which is the combination appearing in quantities computed in full Nf = 2 W%PT. The chiral 



phase diagram for Wilson- type fermions depends on the LEC c 2 (55) and in particular two 
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FIG. 4: Parameter space of the results obtained fitting the numerical data of ref. [7] using eq. (CI ). 
See the main text for details on the fit performed. Only results where a reasonable fit is obtained 
are showed. 



scenarios can take place depending on its sign. In our case only one of the two LECs appears 
in the discretisation correction to the spectral density. For large number of colours N c , one 
finds Wq/Ws ~ l/N c [12], indicating that W$ is the dominant contribution in the coefficient 
c 2 . This implies that independently on the sign of W$ a positive value for W$ would suggest 
a negative value for c 2 . 5 Assuming for instance Wq = W$/3, then the results we obtain 
from the fit A = -0.55 imply that a 2 W$ = 5.4 • 10 6 MeV 4 which for a = 0.08 fm would 
correspond to c 2 — —(473 MeV) 4 . It would be interesting to extend the numerical data in 
order to establish more accurately the sign of and the absolute value of Wg. This will give us 
valuable informations on the behaviour of the spectrum of the Hermitean Wilson operator 
and a hint on the scenario for the chiral phase diagram that takes place with clover fermions 
and Wilson gauge action. 6 

We close this section by stressing again an important remark. Our results seem to indicate 
that apart from the quark mass, the lattice spacing and the lattice size, there is an additional 
scale entering in this problem, namely the eigenvalues of the massless Dirac operator. In 
particular, after the analytic continuation of eq. (|9| the mass ^Jm 2 + /z 2 becomes y/X 2 — m 2 , 
which parametrizes the "distance" from the threshold of the spectral density. Even if we 



5 While for other discretisations existing numerical results of refs. [401 - 142] are consistent with ci < 0, for 
non-perturbative 0(a)-improved Wilson fermions and Wilson gauge action the value of ci is still unknown. 

6 If the negative sign of C2 found in this work will be confirmed by other studies, it will imply that the pion 
mass never vanishes. With our numerical estimate of C2, the minimal pion mass at a = 0.08 fm would be 
around 130 MeV, which is well below the values actually simulated with this action. 
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assume, as it is done in this work, that all masses are in the p— regime the value of this 
additional scale can be arbitrarily small. In order to get meaningful results we have to 
implicitly assume that this parameter obeys the same power counting as the quark masses 
with respect to the lattice spacing a and the lattice size L. While this fact seems trivial 
at this stage, it is not totally obvious to infer it from the partially quenched initial setup, 
where valence quarks are introduced as probes to obtain the spectral density. This implies 
that, at fixed a and L, our results will not be valid in the vicinity of the threshold, since 
from one side the p-regime expansion will fail, and from the other side the GSM counting 
will not be valid anymore, because close to the threshold, as we have discussed before, the 
NLO corrections induced by the finite lattice spacing dominate the LO result. 

We stress that these two effects are decoupled: close to the threshold, there is a scale 
which in continuum at finite volume needs to be treated with the e-expansion, and 
which in W%PT in infinite volume needs to be treated with the Aoki power counting. The 
way the two effects combine close to threshold is a not trivial issue and it deserves further 
investigations. 

We conclude that our formula is a useful tool, combined with numerical data, to un- 
derstand the behaviour of the spectral density of the Hermitean Wilson operator. We also 
consider this work a necessary and important step towards a complete theoretical under- 
standing of the behaviour close to the threshold of the spectral density of the Hermitean 
Wilson operator. 
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Appendix A: SU(m|n) graded group: conventions and properties 

In this Appendix we collect conventions and properties related to SU(m|n) [291 03] which 
are relevant for the partially quenched Chiral Effective Theory in the graded-group formu- 
lation. 



1. Supermatrices 

For n, m positive integers, a square even supermatrix is defined as a (m + n) x (m + n) 
matrix with the structure 

jj ( A m xm B n xm \ ( a i \ 

\ n n r [ } 

\ ^mxn J-^nxn- / 

where A and D have elements in the even subspace of the Grassmann algebra, while B 
and C have elements belonging to the odd subspace of the Grassmann algebra (see [13] 
for definitions and properties of the Grassmann algebras). We define the supertrace of a 
supermatrix 

Str([/) = Tr(A) - Tr(L>), (A2) 
and the superdeterminant of an invertible supermatrix as 

Sdet(ET) = Det(A - BD^^/Det^D). (A3) 



2. Superalbegras 

We now consider a square (m+n) x (m+n) matrix with complex entries and the structure 

M=( ^ mxm ^ nxm ) . (A4) 

M is said to be even if B = C = 0, whereas it is ocM if A = D = 0. The degree or parity of 
M is defined to be 

f if M is even . . _ 

deg M = < . r „ , . , , A5 
6 \ 1 if M is odd y ' 

One can show that set of all complex linear combinations of these matrices form an associa- 
tive superalgebra. The supertrace can be defined analogously as in the case of supermatrices, 

Str(M) ee Tr(A) - Tr(£>). (A6) 



3. Lie supergroups and superalgebras 

The Lie Supergroup SU(m|n) is the set of (m + n) x (m + n) even supermatrices [/ 
satisfying the conditions 

U ] U = 1, Sdet([/) = 1. (A7) 
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By parametrizing U as 
these conditions convert into 



U = e i(p , 



= <t>\ Str(0) = 0. 



(A8) 



(A9) 



The Lie superalgebra associated to SU(m|n) can be constructed by supplementing the asso- 
ciative superalgebra of (m + n) x (m + n) complex matrices defined in the previous subsection 
with the generalised commutators 



[M,N] = MN - (_l)(de g M)(deg7V) ArM? 

and by requiring the conditions 

M = M\ Str(M) = 0. 



(A10) 



(All) 



A basis for the Lie superalgebra of SU(m|n) can be chosen to be a set of Hermitian matrices 
T a = ( Ta y ( a = i ; . . . ) ( m + n )2 _ !) ; a definite parity 



and zero supertrace 

The normalisation is given by 



deg(a) = degT a , 
Str(T) a = 0. 



Str(T a T 6 ) 



9 



ab 



where 



/I 



9 ab = 



\ 



2 ' 



1, . . . , m 2 



(A12) 
(A13) 
(A14) 



m 2 , . . . , m 2 + 2mn — 1 



m 2 + 2mn, . . . ,(m + n) 2 — 1, 



(A15) 



-1/ 



with g ab = (— l) dc s( a ) dc s( fe )g 6a . The elements from m 2 to m 2 + 2mn — 1 are odd matrices, 
while the remaining ones are even. 

In this work we consider the particular case of SU(4|2). With our conventions, T 1 , . . . , T 15 
are the (even) generators of the SU(4) subgroup that includes sea and valence quarks, 
T 32 , . . . ,T 34 (even) are associated with the SU(2) ghost sector, T 16 , . . . ,T 31 (odd) mix the 
ghosts and the quarks, and finally T 35 (even) is a diagonal matrix with components in both 
quark and ghost sectors. 
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For the computation of the propagator it is convenient to introduce the tensors 



k ab = I _}_ g aU + 1 g al5 _ _L g *A ( * W4 + 1 615 _ J_ g b3 5 ) ? (Alg) 



V~r ' Vq* sir ) VVr ' Vq" y/r 
hab = , + + ^ ^ + ^ + ^ 635 , (A17) 



which satisfy 



35 



^ h ac k cb = 0. (A18) 



c=l 



Appendix B: The spectral density in W%PT with the replica method 



In this appendix we present the same calculation done in sec. IV applying the replica 
method [28J in alternative to the graded group method. We consider a partially quenched 
theory with N s sea quarks of mass m s and N v quenched valence quarks of mass m v . For 
simplicity we set the volume to infinity; finite-volume corrections can be computed in a 



standard procedure, like explained in sec. IV C The replica method consists in enlarging 



the valence sector to k replica of N v valence quarks. The full symmetry group at zero 
quark mass is therefore SU(A r s + N r )LxSXJ(N s + N r )n, with N r = kN v , and the generating 
functional at the quark level is given by 

Z PQ ,re P uca(J) = J [dA„]Det(£> +m v + J)^Det(L> +m v ) N ^Det(D + m s ) ^ e ~ s ^ , (Bl) 

where we have introduced source terms J. In the limit N r = one reproduces the generating 
functional in the graded group method j26l EZj 

Z PQ , 9 raUJ) = J W ^p^y' DetP + m.) N '<r 8 ' {A >\ (B2) 

which is obtained in a SU(A S + N v \N v )lxSXJ(N s + N v \N v )r theory, where internal valence 
quark loops are cancelled by introducing quarks with the wrong statistics (ghosts). This 
equivalence is translated into the chiral effective theory and has been verified at the per- 
turbative level [28]. In this appendix we explicitly show how this equivalence is realised for 
the computation of the spectral density starting from the partially quenched pseudoscalar 
density. 

For our specific computation we consider the case N s = 2, with a doublet of twisted 



valence quarks. The chiral Lagrangian is the same as in eq. (16), with the "Supertrace" 
replaced by the conventional trace over the group SU(A S + N r ). The mass matrix has the 
form 

Ai = diag( ^J^s ' 171 + ^/i^r 3 , . . . , m + ifi v r 3 ). (B3) 

2x2 sea Nv=2 

^ v ' 

N r =kN v 
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The pseudo Nambu-Goldstone fields are parametrised by 

U(x) = u v e 2i ^ x)/F u Vl ^ = Y^CT a , (B4) 

a 

where now T a , a = 1, . . . , (2 + N r ) 2 — 1 are the generators of SU(2 + N r ), with the normali- 
sation convention 

Sjab 

Ti(T a T b ) = — , (B5) 

while uy represents the ground state of the theory. By minimising the LO potential we 
obtain 

u v = diag( ^ , , (B6) 

2x2 sea N r 



with the definitions already given in eq. (15) 



sin wo = — — , cos wo = , m P = x/m 2 + (B7) 

mp mp 

In this framework it is convenient to write down the pseudo Nambu-Goldstone propagator 
in explicit components 

(U{x)U{y)) = \ [5c b 5 da G\x - y, M 2 ab ) - S ca 5 db E(x - y, M 2 aa , M 2 CC )} , (B8) 



with a, b, c, d — 1, . . . , (2 + N r ). The infinite-volume scalar propagator is defined in eq. (21 ), 
while 

1 f e ipx 
E( X - y, M L Ml) , ^ j *P v + MLW + Mi)F(p y ^ 

2 N 

f ^-WTm!, + WTM' (mo) 

with 

( M 2 S = 2mB if a, b = 1,2 

M 2 b = B(m a + m b ) = I M 2 V = 2m P B if a, b = 3, . . . , N r (Bll) 

{ M 2 V = {m + m P )B if a{b) = 1, 2; 6(a) = 3, . . . , N r . 

The observable we have to consider in order to extract the spectral density is the expectation 
value of the pseudoscalar density defined in eq. ([8]). The pseudoscalar density can be obtained 
by taking the functional derivative of the action with respect to appropriate sources, which 
are introduced by the following procedure: 

M -)• M + fV(x), M ] -> M ] - fV(x), (B12) 

where f 3 has non-zero elements only in one of the replica of the valence sector 

f 3 = diag(^0^ Jv r^,0). (B13) 

2x2 sea N v =2 

N r 
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At LO we obtain the (continuum) expectation value 

(V 3 ) LO = 2i£sinw - (B14) 
This result is independent on N r , and coincide with the one obtained with the graded group 



method, eq. (28). At 0(p ) we obtain 



{V 3 )nlo 



2iS sin(w ) < 1 + 5 cot tu — 



1 



2G 1 (0,M 2 ,) + A r G 1 (0,M, 



2 ) 

w I 



~ E{0, Ml, Ml) - 16L 6 (2Ml + N r M 2 J - AM 2 UV (H 2 + 2L 8 
- 8a (W e (2 + N r cos u ) + W 8 cos u ) 



(B15) 



The shift 6 is the 0{p 2 ) correction to the ground state angle uq that must be computed by 



minimising the NLO potential. Unlike the graded group case discussed in sec. |IV A[ with 
the replica method the minimisation is trivial, since the fields belong to the conventional 
SU(iV s + N r ) group. By taking only the linear term in N r we get 



16a sin coo J ( „ r W 8 \ t 2a ( Wi 



(B16) 



which coincides with the shift calculated with the graded groups in sec. IV A, eq. (44). The 



final result in the limit N r — > for the partially quenched pseudoscalar density at NLO is 
then 



{V 3 )nlo 



2iS sin wo s 1 + 5 cot uq + 



F 2 



-G\0, Ml) - 2G 1 (0, Ml) + 32L 6 M 2 S (B17) 



1 



(Ml - M s 2 JG 2 (0, Ml) + AMl(H 2 + 2L 8 ) + 8a (2W 6 + W 8 cosu ) 



where we have used 



E(0,Ml,Ml)\ Nr ^ = ±GHQ.Ml) -^iMl MDG^Q. Ml). 



(B18) 



This result for the pseudoscalar density is fully equivalent to the one obtained with the 



graded groups, eq. (29), as expected 



Appendix C: Some explicit formulae 

In this Appendix we collect some explicit formulae which we omitted in the main text 
for clarity reasons. 
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1. Integrated spectral density at NLO 



We report the integrated spectral density computed at NLO in W%PT. Starting from the 



result presented in eq. (46) for the spectral density we obtain: 



M 2 

N(A 1 ,A 2 ,m) NLO = / d\[p Q (X,m) + p Q (-X,m)} NLO 

2S 

71 



Am 2 



(CI) 



+ 



(4tt) 2 F 4 
+ (A 2 , — 2m 2 ) arctan 



- ^(A 2 - A 2 ) + m(l + 3L 6 ) (^Af. - m 2 - ^/a 2 



- - m 2 



(A x — 2m ) arctan 



v/Af 



2 _ m 2 



ni 



m\ A 2 , — m 2 21o, 



+ mWAf - m 2 21o, 




where A is given in eq. (47). 



2. Finite volume corrections 



In this section we give explicit expressions for the finite volume corrections defined in 



eqs. (53), (54). It is convenient to represent the finite-volume propagator g r (M 2 ) [30j defined 



in eq. (20) as sum over modified Bessel functions K v : 

1 



9r(M 2 



r(r)(4vr) 



{nl,n2,n3,n4}^0 



,||,M 2 



(C2) 



with ql = {{n\ + n 2 2 + n\)L 2 + n\T 2 ) and 

F u (a,z) =2 2 K v {2^). 



(C3) 
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Using this representation one obtains the finite-volume correction to the spectral density 



2A 



[Apl(\,m) + ApV(-\ 7 m)} NLO = - j=— £ 



{nl,n2,n3,n4}^0 



Re 



\2F 2 



, iV A 2 — m 2 



- 2Re 
+ VA 2 — m 2 Im 



F_! (||f,m WA 2 -m 2 



+ mRe 



(C4) 



The corresponding correction for the integrated spectral density N V (A 1 , A 2 , m) NLO for Ai = 
m is given by 

M2 OV2 /A2 _ -_2 

A^(m,A 2 ,m) rao = y [Apg(A, m) + Ap$(-A, m)] NLO d\ = — V * (C5) 



X! i /T2 

{nl,n2,n3,n4}^0 I V Ji 2 



zlm 



F - 2 (§W A2 - m2+m 



V4F 2 



m 



Jm 



V 2F 2 



- m 2 



+ Re 



F-i(^,<A/A2-m=» 
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